Charge-fluctuation contribution to the Raman response in superconducting cuprates 
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We calculate the Raman response contribution due to collective modes, finding a strong depen- 
dence on the photon polarizations and on the characteristic wavevectors of the modes. We compare 
our results with recent Raman spectroscopy experiments in underdoped cuprates, La2_xSrxCu04 
and (Yi,97Cao.3)Ba2Cu06.05, where anomalous low-energy peaks are observed, which soften upon 
lowering the temperature. We show that the specific dependence on doping and on photon polar- 
izations of these peaks is only compatible with charge collective excitations at finite wavelength. 
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There is an increasing experimental evidence that in 
the high Tc superconducting cuprates there are finite- 
energy excitations, which are distinct from the usual 
single-particle excitations characterizing the spectra of 
standard metals. In optical conductivity there are 
many examples of absorption peaks at finite frequen- 
cies, which are quite distinct from the zero-energy Drude 
peak, occurring in hole-doped [La2-xSrxCu04 (LSCO) 
m, Lai.6-vNdo.4SrxCu04 (LNSCO) 0, BiaSraCuOe 
(Bi2201) |3] as well as in electron-doped materials 
[Nd2_xCexCu04 (NCCO) Hi]. In aU these cases absorp- 
tion peaks are observed below 200 cm~^ which substan- 
tially soften upon decreasing the temperature T . It is 
natural to assign to these excitations a collective char- 
acter. Moreover in the case of Bi2201 there is a clear 
evidence of scaling properties for the finite-frequency ab- 
sorption 0. This suggests that these excitations are 
nearly critical fluctuations associated to some critical- 
ity occurring in the cuprates below a critical line T*{x), 
ending near optimal doping at T = 0. The proposal of 
a quantum critical point in the str ong ly underdoped 
or near optimal doping 0, B S ll^ llH is by now acquir- 
ing consensus in the community, but the nature of the 
ordered phase still remains debated. In this regard it is 
quite important to establish whether the collective modes 
(CM) responsible for the finite-frequency absorptions oc- 
cur at finite wavevectors and are therefore associated to 
some form of spin and/or charge spatial ordering 0,0 
or have infinite wavelength |3,|9|- 

In order to clarify this issue, in this letter we calcu- 
late the Raman response contribution due to collective 
charge-order (CO) excitations associated to stripe fluctu- 
ations and eventually use our results to interpret recent 
findings of Raman-scattering experiments 0, 0| . Ra- 
man scattering is an important tool in the above issue for 
several reasons. First of all it is a probe directly access- 
ing the bulk properties in the relevant frequency range. 
Most importantly, a suitable choice of the polarization of 
the incoming and outgoing photons selects the regions in 
momentum space from where excitations are originated 
[l^ . As a consequence, characteristic features arise in 



the Raman response with distinct doping and polariza- 
tion dependencies (12„,13]. which in turn allow to extract 
valuable informations about the momentum dependence 
of the excitations. This is the main issue of this work. 
We consider the doping-dependent fine Tco{.x) ^ T*{x) 
[isj . In the underdoped regime below this line the system 
would order in the absence of competing effects. There- 
fore strong CO fluctuations are present in a large portion 
of the phase diagram and provide an additional chan- 
nel for the Raman response, besides the usual one ob- 
tained from Fermi- liquid quasiparticles (QP). The role 
of fluctuations in the response functions was established 
by Aslamazov and Larkin for the paraconductivity in 
superconductors This theory was extended to the 

particle-hole channel for one-dimensional charge-density- 
wave systems in Refs|l3[li In this latter case a partial 
cancellation leads to a less singular paraconductivity as 
compared to the particle-particle case. Here we recon- 
sider this scheme in two dimensions for the Raman re- 
sponse diagrammatically represented in Figs. 1(a) and 
1(b) involving the incommensurate charge CM. We find 
that the above cancellation does not generically occur 
for the corresponding Raman diagrams, generating im- 
portant contributions, which in turn depend on the se- 
lected symmetry. We will argue that only these contri- 
butions account for the strong symmetry dependence of 
the Raman response. Other diagrams with self-energy 
and vertex corrections are weakly anisotropic and pro- 
vide the usual weakly temperature-dependent QP con- 
tribution. Specifically each dashed line appearing in the 
Raman response in Fig. 1 (a,b) is the gaussian diffusive 
CM propagator in Matsubara frequencies 



I?(q,[j„i) = (|cj,„| + riq)" 



(1) 



where fiq — v\(\ — qd + rn(a;, T) with v a constant elec- 
tronic energy scale (we consider a unit lattice spacing). 
m{x, T) is the mass of the CM and encodes the distance 
from the critical line Tco{x). This propagator is domi- 
nant at zero frequency and q = Qc, the wavevector set- 
ting the modulation of the most singular charge fluctua- 
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FIG. 1; Direct (a) and crossed (b) diagrams for the fluctu- 
ation contributions to Raman spectra. The dots represent 
Raman vertices. Sohd lines represent fermionic QP propaga- 
tors and dashed lines represent CM propagators, (c) and (d): 
Hot spots joined by the critical wavevector qc = (0, —Qc) rep- 
resented by arrows. Equivalent hot spots have arrows of the 
same type (solid or dashed). The -I- (— ) sign mark the regions 
where the Big and B29 vertices are positive (negative). 



tions. The fermionic loop in Fig. l(a,b) has the form 

7a,/3(k)G(k,e„ + fli) 

n k 

X G(k - q,£„ - Lj„i)G(k,e„), (2) 

where 7a,/3(k) = d^E\^/ dkadkp and C is a constant de- 
termined by the coupling of the Raman vertex with the 
incoming and outgoing photons and the coupling g of 
the CM with the fermions. The choice of the a and (3 
components depends on the polarization of the incoming 
and outgoing photons A suitable choice of these 

polarizations corresponds to specific projections of the 
7(k) vertex on cubic harmonics of the square lattice. 
To start our discussion we chose the polarization corre- 
sponding to a Raman vertex with Big symmetry = 
coskx — cosfcy vanishing along the (0,0) (±7r, ±7r) 
directions. We are interested in the dominant contribu- 
tions of the diagrams in Fig. l(a,b), which occur when 
the CM are around qc. Therefore we set q = qc in 
Eq. In the k summation the largest contribution 

is obtained when the three fermions are around the so- 
called hot spots, that is those regions on the Fermi sur- 
face which can be connected by qc. For each given qc in 
Eq. (I2Jl, to avoid cancellations, the sum over k must en- 
counter equivalent hot spots where the vertex jBigik) 
does not change sign. Since jBig changes sign under 
the zbcc vs. ±?/ interchange, this condition is fulfilled by 
stripes or eggb ox CO fluctuations at not too small dop- 
ing HEEHmi, where qc ~ 27r(±0.2, 0), 27r(0, ±0.2) 
[see Fig. 1(c)]. On the other hand, when for the same 
qc the i?2g vertex 75^ — sin k^ sin ky is considered in 
Eq. ^ , the leading contribution from hot spots vanishes 
since the equivalent "available" hot spots [see Fig. (Id)] 
give contributions from regions where jB2g is opposite 
in sign. With similar arguments one can show both in 
the superconducting and in the normal phase, that 
spin CM at qc = (±7r, ±7r), although having the simi- 
lar hot spots as the CO ones, give rise to non-vanishing 



vertices for the Aig symmetry only. The same holds for 
spin CM with diagonal incommmensurate wavevectors 
qc ~ (±(7r — (5), ±(7r — SJ), which is the case of LSCO at 
low doping (x < 0.04) j23. On the other hand, at higher 
doping, for small vertical and horizontal incommensura- 
tion qc = (±(7r — (5), ±7r), (±7r, ±(7r — 6)), there is no 
complete cancellation in the Big symmetry. This small 
contribution could add to the CO contribution that we 
are going to evaluate from the diagrams of Fig. l(a,b), 
but cannot explain the anomalous Raman response in the 
deeply underdoped regime. 

The resulting CM contribution to B 
man scattering is given by 
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where h{z) is the Bose function. 
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and z± = {z ± w/2)(l + {z ± w/2)Vt^o)- Here ~ 
100 — 500 cm~^ is an ultraviolet cutoff of the order of the 
frequency of the phonons most strongly coupled to the 
electrons and driving the systems near the CO instabil- 
ity 0] . The above expression of Ax" is actually calcu- 
lated by considering CO collective modes with a semiphe- 
nomenological spectral density of the form 
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In the limit of infinite wq, one recovers the spectral den- 
sity of the above critical diffusive propagator Z3(q, Wm). 

Experimentally [l^ in Lai.9Sro.iCu04 an anomalous 
Raman absorption is observed in the Big symmetry, 
while the behavior of the B2g spectra can be accounted 
for by QP only. Moreover the anomaly in Big softens 
upon lowering temperature and tends to saturate sug- 
gesting a quasi-critical behavior. In Fig. 2(a) we report 
a comparison between our theoretical calculations and 
experimental data. According to the neutron-scattering 
experiments in this material jlil l2l| we choose qc along 
the (1,0) or (0,1) directions, giving a non- vanishing vertex 
for the Big symmetry only. We adjust at one tempera- 
ture the overall intensity by choosing the vertex strength 
Ks^g of Eq. lO and the ultraviolet cutoff wq to reproduce 
the lineshape. Then, keeping fixed these parameters at 
all temperatures, we tune the mass m{T) to reproduce 
all the curves. The agreement is manifestly satisfactory, 
indicating that the strong temperature dependence of the 
Raman absorption peaks is basically ruled by the tem- 
perature dependence of the low-energy scale m{T) and by 
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FIG. 2: (color online) (a) Comparison between experimen- 
tal and calculated Raman spectra on Lai.goSro.ioCu04. The 
intensity is chosen to reproduce the data with A^^^^ = 1.7, 
the high-frequency cutoff is ujo = 250 cm~^, and the mass is 
reported in the inset, (b) Same as (a) for Lai.98Sro.o2Cu04 
with A|, =0.85. 



the temperature dependence of the Bose factors. The in- 
set reports the CM mass m(x = 0.10, T) needed to fit the 
experimental curves: At large to moderate temperatures 
m{T) has a linear part, which has an intercept on the T 
axis at a temperature T* of the order of 70 K. At lower 
temperatures m tends to saturate and the system crosses 
over to a nearly ordered regime with a finite m{T). This 
behavior is clearly consistent with the behavior expected 
for the mass of critical modes in the underdoped region: 
At large temperatures the system is in the quantum crit- 
ical regime with m oc {T — T*) , while below a crossover 
temperature T* [which for LSCO at a: = 0.10 is indeed 
about 60 - 80K, see, e.g.. Fig. (la) of RefEi] the CO 
transition is quenched by other effects. 

If the CM had wavevectors along the diagonal (1,1) and 
(1,-1) directions and small modulus, the role of the Big 
and B2g would be interchanged, with this latter display- 
ing the anomalous Raman ab sorp tion. Noticeably, exper- 
iments in Lai.98Sro.o2Cu04, [l^l, where neutron scatter- 
ing detect spin incommensurate (but likely also stripe) 



order along the diagonal directions '23!|, do show that 
the anomalous Raman absorption is present in the i?2g 
symmetry and is absent in the Big. Fig. 2(b) reports a 
comparison between experimental data of Ref. [12] and 
a calculated CM contribution for this B2g case. In this 
case, consistenly with the expected high value of T*, the 
temperature dependence of the mass no longer displays 
a clear linear behavior and therefore the estrapolation of 
the high-temperature "linear" part to identify T* is not 
well justified. Nevertheless, as a crude estimate from the 
last three points we obtain T* of the order of 150 K, 
which is again consistent with the T* values reported in 
the literature (see, e.g., Ref. Iisi) . Notice also that the 
tail of the data at the lowest temperature (T = 88 K), 
substantially lower than T*, is not well reproduced by 
our calculations based on a diffusive form of the CM. 

A similar behavior of the Raman spectra is observed in 
(Yi.g7Cao.3)Ba2Cu06.o5 sample 13] at filling p ~ 0.03, 
with a temperature-dependent peak in the i?2g symmetry 
only. Again it is quite natural to attribute this peak to 
stripe fluctuations along the (1,±1) directions and long 
wavelength [qc ~ (p,p)]. Fig. 3 reports the comparison 
between experiments and our theoretical analysis. Again 
a value of T* ^ 120 K may be estrapolated from the high 
temperature mass behavior. 

The data for YCBCO also display a remarkable change 
in the lineshape at the lowest temperatures (more pro- 
nounced than in Lai.98Sro.o2Cu04 at T — 88K. In par- 
ticular the peak at T = 11.0 K is much narrower and 
decreases more rapidly at high frequencies. This change 
in the lineshape naturally marks a change in the nature 
of the charge fluctuations from the the high-temperature 
diffusive behavior in the quantum critical regime to a 
sharper mode behavior typical of the low-temperature 
underdoped regime. In this case the diffusive form of 
the propagator fails in reproducing the lineshape, while 
a good description may be obtained by a factorized form, 
according to previous treatments of the nearly ordered 
underdoped regime 0|. Therefore we also carried out 
the calculations using 



i)(q,CJ„) = W{iujn)Jxiq)Jyiq), 



(6) 



where Jx,y{<l) = J^l [7^ + - fc,y)^] • Here M is 
a normalization factor and the inverse correlation length 
is 7. The (real-frequency)-dependent part W{uj) = 
/ dvA{v)2v/{uj'^ - z/2) with A{u) ~ r/[(a; - wcm)^ + F^] 
is a normalized lorentzian centered around the typi- 
cal CM energy ojcm with halfwidth F. now 
reads Axg = B{T)J[uj/V,T/T,ujcm/T) with B{T) oc 
9{TY/b{Tf^)], 9 being the coupling between the CM 
and the QP [25j. As shown by the dashed curves in Fig. 
3, the factorized "propagator" provides a good descrip- 
tion at the lowest temperature, while in the intermediate- 
temperature regime (T = 86.2 T = 127. iK) the 
diffusive propagator becomes progressively a better de- 
scription of the charge fluctuations. We also notice that 
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FIG. 3: (Color online) Comparison between experimental and 
calculated Raman spectra on (Yi.97Cao.3)Ba2Cu06.05 with 
diffusive propagators [Eq. (solid lines). The intensity is 
chosen to fit the data with A^^^ = 0.41, the high-frequency 
cutoff is LOO ~ 100 cm~^, and the mass is reported in the in- 
set. The dashed lines are obtained with factorized "propaga- 
tors" [Eq. ©] with overaU intensity B{T) = 28, 27, 19, 12 
(cps/mW) for T = 11 K, 55.2 K, 86.2 K, 127.3 K. r(T) 
and LOOM are also reported in the inset. 



the mass m(T) of the quantum-critical diffusive prop- 
agator smoothly connects with the energy of the low- 
temperature modes: m{T) ^ ujcm- 

As far as LSCO samples are concerned, the diffusive 
form of the modes seems to persist down to the lowest 
temperatures in the available published data. We notice 
in passing that the two values for T* obtained at x — 0.02 
and x = 0.10 linearly estrapolate to T* = at a; = ~ 
0.17, consistent with the theoretical position of the CO 
quantum critical point |^ and not inconsistent with 
most of the experiments on T* in LSCO 11]. 

In order to emerge clearly from the other electronic 
excitations, the CM's need to have rather small masses 
(from the insets of Figs. 2 and 3 a few tens of cm~^). 
This is the case in underdoped systems where the fi- 
nite and sizable values of T* carry to finite tempera- 
ture the small values of the mass, which should vanish at 
Tco{x) ^ T* if a genuine long-range order would occur. 
Around optimal doping, where Tcoi^c) = 0, anomalous 
peaks are hardly visible because superconductivity inter- 
venes at higher temperatures before the mass, being pro- 
portional to T, becomes sufficiently small. Similarly we 
argue that the observation of CM peaks in Bi-2212 sys- 
tems is unlikely. In these materials, STM experiments 
suggest that the CO coherence lenght is smaller than 5-6 
lattice spacings indicating large mass values, in con- 
trast to underdoped LSCO materials, where according 
to neutron scattering experiments the strijse correlations 
may extend over tens of lattice spacings 19 . 

In summary we calculated the CO-CM contribution 
to the Raman response function. The idea that charge 
modes are important in the low-energy Raman spectra 



is also supported by recent experiments in ladder com- 
pounds |2g. Here we showed that the CM contribu- 
tion quantitatively accounts for the observed anomalous 
peaks in Raman spectra of LSCO and YCBCO. In this 
explanation a crucial role is played by the symmetry of 
the Raman vertices: We demonstrate the stringent con- 
nection between the symmetry properties of the Raman 
spectra and the doping, which implies a dependence of 
the CM on finite wavevectors specific of stripe or eggbox 
fluctuations as inferred from neutron scattering. This 
rules out for the interpretation of these Raman exper- 
iments local (i.e., at all momenta) excitations like po- 
larons, disorder-localized single particles or resonating- 
valence-bond singlets. Also spin CM do not comply 
with these symmetry requirements over the whole doping 
range examined above. For the same reason other excita- 
tions peaked at q = (like superconducting pair fluctua- 
tions or time-reversal-breaking current fluctuations) are 
not appropriate. Although we cannot a priori exclude 
that the proposed CO collective excitations are concomi- 
tant with these other more elusive forms of criticality, at 
the moment the anomalous Raman absorption can only 
be interpreted with CM peaked at finite momenta. 
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